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In a recent study [3| the stationary state of a parallel-update TASEP with varying system length, 
which can be regarded as a queueing process with excluded- volume effect {exclusive queueing process, 
EQP), was obtained. We analyze the dynamical properties of the number of particles {Nt) and 
the position of the last particle (the system length) (Lt), using an analytical method (generating 
function technique) as well as a phenomenological description based on domain wall dynamics and 
Monte Carlo simulations. The system exhibits two phases corresponding to linear convergence or 
divergence of (Nt) and (Lt). These phases can both further be subdivided into high-density and 
maximal-current subphases. The predictions of the domain wall theory are found to be in very good 
agreement quantitively with results from Monte Carlo simulations in the convergent phase. On the 
other hand, in the divergent phase, only the prediction for (Nt) agrees with simulations. 

PACS numbers: 02.50.-r, 05.70.Ln 



I. INTRODUCTION 

Queueing processes have been studied extensively, es- 
pecially due to their practical relevance [2^4}]. However, 
usually the spatial structure of the queues is neglected 
and the particles in the queues do not interact with each 
other. On the other hand, the totally asymmetric sim- 
ple exclusion process (TASEP) which has a spatial struc- 
ture and excluded-volume effect (hard-core repulsion) is 
one of the best-studied interacting particle systems _5]. 
Nowadays the TASEP is a basic model for pedestrian and 
traffic flows [i,0|. 

Recently a queueing process with the excluded-volume 
effect, the exclusive queueuing process (EQP), was intro- 
duced in P and @ independently, where the model was 
formulated as continuous-time and discrete-time Markov 
processes, respectively. This model can be rephrased as 
the TASEP on a semi-infinite chain with a new boundary 
condition, see Fig. [TJ The left end is interpreted as the 
end of the queue where new customers arrive. Therefore 
particles can enter at the left site next to the leftmost 
occupied site. The (fixed) rightmost site corresponds to 
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FIG. 1: Exclusive queueing process. 



the server. Here particles can leave the system after get- 
ting service. In the bulk a particle can hop to its right 
nearest neighbor site, if the target site is empty. After 
the study |9|, where the bulk hopping rule is determinis- 
tic, the model with discrete time and probabilistic bulk 
hopping was analyzed [l| . 

Earlier works [ll, H, focussed on the exact probability 
distribution, physical quantities in the stationary state, 
and the conditions under which the EQP converges to the 
stationary state. In this paper, we study the dynamical 
properties by considering the number of particles and 
the system length which is defined as the position of the 
leftmost particle. We use the same formalism as in [H, |1] , 
i.e. discrete time and parallel-update scheme. For generic 
hopping probability p, we will introduce a domain wall 
prediction, checking it by Monte Carlo simulations. In 
the deterministic hopping case p = 1, a, rigorous analysis 
is available, by using the generating function technique 

This paper is organized as follows. In SecHH we define 
the model as a discrete-time Markov process, and review 
its stationary state based on Jj , which can be generalized 
to the inhomogeneous injection case, see also App. [X] In 
Sec. HID we introduce a phenomenological argument on 
how the number of particles (Nt) and the system length 
(Lt) converge or diverge, showing simulation results. In 
Sec. II V[ we derive the asymptotic behaviors of {Nt) and 
{Lt) rigorously for p = I, imposing the initial condition 
that there is no particle on the chain. Section |V] is de- 
voted to the conclusion of this paper. 



II. MODEL 
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The EQP is defined on a semi-infinite chain where sites 
are labeled by natural numbers from right to left (Fig.[T|). 
Particles can enter the chain with probability a only at 
the left site next to the leftmost occupied site. A particle 
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hops to its right nearest neighbor site with probabihty 
p, if it is empty, and exits at the right end of the chain 
with probabihty /3. If there is no particle on the chain, a 
particle enters at site 1 with probability a. These transi- 
tions occur simultaneously within one time step, i.e. we 
apply the fuUy-parallel-update scheme. 

We formulate the EQP as a discrete-time Markov pro- 
cess on the state space 

S = {0, 1, 10, 11, 100, 101, 110, 111, 1000, . . .} (1) 

where and 1 correspond to unoccupied and occupied 
sites, respectively. In particular, denotes the state in 
which there is no particle on the chain. To simplify the 
notation we do not write the infinite number of O's lo- 
cated left to the leftmost 1. 

Let us review the matrix product stationary state [1] of 
the model, which is a simple extension of that for systems 
with a fixed system length [10|. When 




for /3 > 1 
for ^ < 1 

the stationary state can be expressed as 

p(0)-4, 
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(4) 

Xi ~ D and Xq = E are matrices, {W\ is a row vec- 
tor and \V) is a column vector satisfying the algebraic 
relations 

EDEE ^{l- p)EDE + EEE + pEE, 

EDED ^EDD + EED + pED, 

DDEE =(1 - p)DDE + (1 - p)DEE 
+ pil-p)DE, 

DDED =DDD + (1 - p)DED + pDD, 
DDE\V) =(1 - P)DD\V) + {l-p)DE\V) 

+ p{l-P)D\V), 
EDE\V) ={l-l3)ED\V)+EE\V)+pE\V), 
{W\DEE ^{l~p){W\DE, 
{W\DED ={W\DD +p{W\D, 



(5) 



DD\V) 
{W\EE 



_p{l 



j^D\V), ED\V)^^E\V), 



=0, {W\ED =p{W\D, {W\D\V) 



These relations are closely related to those for the sta- 
tionary state of the parallel-update TASEP with ordi- 
nary open boundary condition jlli] . The normalization 
constant is expressed as 



Z = 1 



E 

L>1 

2(1 



p{l - a) 
-a)/3 



R-p + 2{l~ a)l3 



{W\D{D + E)^-^\V) (6) 



(7) 



with R = ^Jp{p~Aa{\ - a)). The average number of 
particles {N) and the average system length (the position 
of the leftmost particle) (L) are calculated as 
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in the stationary state 
on the critical line a = 



Note that 



(N) and (L) 
/? > 1- 



(8) 
(9) 
diverge 



^ yr^), /? > 1 - 

where the stationary state exists. 

A generalization of the model, where the entry prob- 
ability depends on the system length, also has a matrix 
product stationary state, see App. \K\ 
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FIG. 2: Phase diagram for the EQP. The parameter space 
is divided into two regions with and without the stationary 
state. 



III. DOMAIN WALL PICTURE AND MONTE 
CARLO SIMULATIONS 

In this section, we discuss the time evolution of the 
average number of particles {Nt) and the average system 
length {Lt) corresponding to the position of the leftmost 
particle. 

In the ordinary open boundary case, where the length 
of the system is fixed, a domain wall moves rightward 
or leftward, or exhibits a random walk depending on the 
boundary parameters [l^ . In the same way, we will dis- 
cuss how the system length {Lt) moves. We also observe 
how the average number of particles {Nt) changes as well. 
The continuity equation 



{Nt+i) ~ {Nt) = Jl' 



J? 



(10) 



holds, where J™ and J°" are the flows of particles en- 
tering and leaving the system, respectively. The inflow 
J" is always a, which is due to the fact that the site 
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where particles enter is by definition never blocked. In 
other words, our model is not a call-loss system. Under 
the assumption that the outflow J°"* is independent of 
t, we have (Nt) = {a- J°''^)t + (iVo). In fact our simula- 
tions show that both (Nt) and (Lt) decrease or increase 
linearly in time t according to a < ac oi a > ac, respec- 
tively. 



{L,) 



A. Convergent Phase 

When a < ac, the system converges to the station- 
ary state (HI)- We impose the initial condition that 
particles are distributed uniformly with density 



i for /? > 1 - yr^, 

for/3<l-yr^. 



(11) 



As in Fig. [21 (Nf) and (Lt) decrease linearly in time as 

(12) 



a - J°"* 



{Lt 



P 



-t+{Lo), 



(13) 



with 



jout ^ 



1 - v/1 - 4pp(l - p) 



for/3>l-VT^, 
for/3<l-VT^. 



(14) 



Since J™' = 



we have a — J°"* < which means that 



the domain wall moves rightward. 

According to the forms for J°"' and p, we call the 
phases 



MC-C: a < ^ and /3 > 1 - ^1^, (15) 



HD-C: a < 



2 

/?(p-/3) 



and /3 < 1 - Vl -P, (16) 



maximal-current-convergent (MC-C) and high-density- 
convergent (HD-C) phases, respectively, see Fig. 21 

It should be noted that the outflow is given by 
Eqn. (dH) only while < t < J^jL ■ As i ^ 00, the 
outflow approaches a, assuring that {Nt) approaches the 
stationary value ^ . 

Figure [5l shows density profiles in the HD-C and MC-C 
phases. We can observe that the bulk density keeps its 
initial value ([TTl). 



B. Divergent Phase 

When a > Qfc, it is natural to expect that the domain 
wall moves leftward, and the time evolutions of {Nt) and 
{Lt) are expressed by Eqs. and p^ . respectively. 
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FIG. 3: Dynamics in the HD-C (top) and MC-C (bottom) 
phases. The parameter values are chosen as {a,l3,p) — 
(0.2,0.4,0.84) and (0.2,0.8,0.84), and the initial conditions 
as ((A'^o), {Lo)) = (400p,400) and (600p, 600) (with p defined 
by respectively. The green and red lines are data for 

(A'^t) and (Lt) obtained from Monte Carlo simulations, where 
5000 samples are averaged. The black and blue hues corre- 
spond to the predictions of the domain wall theory. Note that 
the asymptotic values are small but non-zero (see Eqs. ([8]) 
and ©): {{N^),{L^)) = (1.70,2.28) and (0.42,0.56), re- 
spectively. 



with the density pT|) and the outflow (|14p . The simula- 
tions imply that this is true for {Nt) with 

{Nt) ^{a-r-')t+{No), (17) 

but fails for (Lt), see Fig. [51 This failure is not unex- 



pected since the predicted velocity 



can be greater 



than 1 whereas the length {Lt) can not be larger than 
t + (Lo) by the definition of the model. However, the 
simulation results (see Fig. [5]) indicate that 



{Lt)^Vt + {Lo), 



(18) 



so that the prediction is qualitatively correct. The veloc- 
ity V satisfies 



V (a -> ac), 
V = l (a = l). 



(19) 
(20) 
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FIG. 4: Subphases of the EQP. 
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FIG. 5: Density profiles {{rj^t) of tlie jth site at time t) in 
the HD-C (top) and MC-C (bottom) phases. The parameters 
(a,/3,p) and the initial conditions {{No), (Lq)) are set to the 
same values as in Fig. O The colored snapshots are obtained 
by averaging 5000 samples of Monte Carlo simulations. The 
black lines represent the predicted densities p according to 
Eqn. (fTTj) . 



and we have exactly (Lt) — t + (Lq) when a = 1. More- 
over, when p — 1, we will show in the next section that 



a~J° 



1 

1+/3 



(21) 



Equations (IT71) and (IT51) can be regarded as the asymp- 
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FIG. 6: Dynamics in the HD-D (top) and MC-D (bottom) 
phases. The parameter values are chosen as {a,j3,p) = 
(0.75, 0.4, 0.84) and (0.75, 0.8, 0.84), respectively, and the ini- 
tial conditions as ((iVo), (Lo)) = (200p, 200) (with p defined 
by HI}). The green and red lines are data for {Nt) and (Lt) 
obtained from Monte Carlo simulations, where 1000 samples 
are averaged. The black and blue lines correspond to the 
predictions of the domain wall theory. 



totic behaviors 



(Nt) 
(Lt) 



(a- 
Vt 



J™*)t + o(t), 
^o(t). 



(22) 
(23) 



In the same way as in the convergent phase, we call 
the subphases 



MC-D: 
HD-D: 



a > 



a > 



2 

P{P - P) 



and 
and 



/3 > 1 - v/T^, (24) 



(25) 



maximal-current-divergent (MC-D) and high-density- 
divergent (HD-D) phases, respectively. Note that the 
densities in the MC-D and HD-D phases are higher than 
(or equal to) ^ and ^3^, respectively. 

It is difficult to predict how {Nt) or {Lt) behaves just 
on the critical line a = ac- For p = 1, however, we will 
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FIG. 7: Phase diagram for p = 1. 



find in tlie next section diffusive behavior on the critical 
line as 



DLVt + o{Vi) 



(26) 
(27) 



with constants Dn and Dl- 



IV. ASYMPTOTIC BEHAVIORS FOR p = 1 

In this section we investigate the asymptotic behaviors 
of {Nt) and (Lt) rigorously for p = 1. Thanks to the de- 
terministic particle hopping we can obtain the generating 
functions of (Nt) and (Lt). For simplicity, we impose the 
initial condition (there is no particle in the system at 
time t — 0), and reset the state space as 



S := {0, 1, 10, 11, 101, 110, 111, 1010, 1011, . . .} 
— {t £ S\t does not contain sequence 00} . 



(28) 



Note that, for p = 1, the sequence 00 never appears if 
the system starts from the initial condition 0. In this 
case, the MC-D and MC-C phases vanish from the phase 
diagram, see Fig. [T] 

We first consider the number of particles, borrowing 
the classification from §| as 



F/'(iV) =Prob 
(iV) =Prob 



^ of particles is N at time t 
A site 1 is occupied at time t 

# of particles is N at time t 
A site 1 is empty at time t 



(29) 
(30) 



for N G Z>o with P/^(0) = 0. These probabilities are 
governed by the following master equation, which was 
found in Q: 



(31) 
(32) 



P,|i(l) =(1 - a)(l - /3)P,^(1) + aPf (0) 

+ (l-a)P,^(l), 
P,-ti W =«(1 - m^'iN - 1) + aPfiN - 1) 

+ (1 - a)(l - /3)F,^(iV) + (1 - a)Pf{N), 
P,l,iO) =(1 - a)P,^(0) + (1 - a)/3P,-^(l), (33) 
P,l,iN) ^a(3P,^iN) + (1 - a)/3P,-^(7V + 1). (34) 

This simplification is due to the deterministic hopping 
p = 1. We choose the initial condition such that 

P„S(0) = 1, p^{N)^P^^{N)^0 {NeN). (35) 



We will check that the average number of particles {Nt) 
converges to the stationary value when a < , and 
show that (Nt) behaves as Eqn. ([^ when a > We 
will also show that (Nt) exhibits diffusive behavior on 
the critical line a 



J_ 

1+/3- 



We define the generating functions of P/^{N) and 
Pf{N) as 



Gf{N) = Y,Pt^{N)z\ 

t>o 

GfiN) = ^P,^(A^)z*, 



(36) 
(37) 



for \z\ < 1. Noting the initial condition ([55]) . we find 
G^(l) - (1 - a)(l - p)zGf{l) + azGf (0) 



+ (l-a)zGf(l), 
Gf{N) = a{l - l3)zG^{N - 1) + azGf{N - 1) 
+ (1 - a)(l - l3)zGf{N) + (1 - a)zGfiN), 



(38) 
(39) 



Gf (0) - 1 = (1 - a)zGf (0) + (1 - a)l3zGf{l), (40) 
Gf (TV) = al3zG^{N) + (1 - a)l3zG^{N + 1). (41) 

From Eqn. PO)) . we have 

l + {l-a)f3zGf{l) 



1 - (1 - a)z 



Inserting this and Eqn. (|41l) into Eqs. 
get a recurrence formula for G^{N) as 



(42) 
and we 



6 



^.A(^. ^ ^ , 1 - (1 - a)(2 - + (1 - a)(l - a - /3 - aP)z^ + (1 - a)^aPz^ . 

^ (1 - a)2/3(l - (1 - a)z)z (1 - a)^/3{l - (1 - a)z)z^ ^ ^ 



(43) 

Gf{N + 1) = + '""'' G^(jV - 1) 



1 - (1 - - /?)z - 2(1 - a)apz^ ^ , , a(l - /?) + a^/^z^^^ 
(l-a)2/3z2 

(=:xGf(iV)+yG^(7V-l)) (TV G Z>2). 



r 



The recurrence formula (|44p has the foUowing solution: are solutions to — xX + y with a; and y as defined in 

Eqn. (gH). Restricting the "initial condition" G^(l) and 



G-(iV) (G-(2) - A_G-(1)) ^^(2) ^"^^ ^^^^ 



A+ - A_ 
where 



\N-1 

- (G^(2)-A+G^(l)) 



(45) 



\Gf(N)\<J2\z\'^j^ (48) 



t>o 



we have 



^ _ 1 - (1 - a)(l - /3)z ~ 2(1 - a)a/3z^ ± r 
^ 2(l-a)2/3z2 G^(2)-A_G^(1) = 0. (49) 



with 



(Note that < |A_| < 1 < |A+|.) From this constraint 
= V(l - (1 - a){l - /?)z)2 - 4(1 - Q;)a/?z2 (47) and the relation G^(l) is determined, and we find 

I 



G^{N) = X^-^G-Hl) = A^-^ + ' ' ' ' "^i; ' " ~ ^~ + (!-(!- ^)^)^ (^.^N), 

2[1 — ay p^{l — z)z^ 

(50) 

Gf (0) = Gf(7V) = (a + (l-«)A_)/3zA^-iG^(l) (iV e N). (51) 

2(1 — Q;)pz(l — z) 

Then we obtain the generating function Gz{N) of the probability that the number of particles is TV as 

G.(0) = Gf(0), G.(iV) = G^(iV) + Gf(7V)-.9(z)A^, (52) 
P-l+{l-a^p + P^~ a/?2) z + (1 - a)(l - /?)/3z2 + ^(1 _ (1 _ z)/3) 
2/3(1 — a)(l — (1 — za)p)(l — z)z 



I 

We also introduce the generating function Q^q of the gen- and Kz of the average number of particles as 
erating function Gz{N) as 

t>o t>o N>a 

z(2a - 1 + z(l - a){l - l3) + r) 
C=i 2(l-z)2(l + z/3) • 



g.C = E Gn{z)C = Go(z) + 5] 5(^)(A-C)^ (55) 

= Go{z) + g{z) ^'^^ The asymptotic behavior of {Nt) is determined by the 

1 ~ A_C degree of the singularity z = 1 of Kz [l3|- When a < 
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1+/3 



we find 



(1 - z)K, 

and (Nt) converges as 

a(l — a) 



a(l — a) 

l3 — a — al3 ' 



(56) 



{Nt 



it 



oo 



(3 — a — afi 

Of course tliis limit value agrees with the stationary value 



with p = 1. When a > we find 



(1 - z)'K, 
and (Nt) behaves as 



1 + 13' 



When a = j^, 



1 + /? 
{a - J°"*)t + o{t) 

we find 



t + o{t) 



{t -> oo). 



and {Nt) behaves as 
(Nt) = 2, 



(58) 

(59) 
(60) 

(61) 



/3t 



^(l + /3)3 



+ o (Vt^ {t oo). (62) 



Now we turn to the behavior of the length of the system 
(the position of the leftmost particle). Let Qt{L) be the 
probability that the system length is L at time t. The 
probability Qt{L) is governed by 

Qt+m - (l-a)Qt(0)+/3(l-a)Qt(l), (63) 
Qt+i{L) = aQt{L - I) + (I - a){l ~ P)Qt{L) 

+ {l-a)PQt{L + l) (64) 

for i G N. The first equation means that, if there is no 
particle at time t + 1, there is no particle at time t and 
no particle enters (with probability 1 — a), or there is 
only one particle on the rightmost site at time t which 
leaves the system and no particle enters (with probability 
(1 — a)/3). The second equation is derived in App. [B] 

In the same way as for the number of parti- 
cles, we define the generating function Mz{L) = 
J2t>Q ^*QtW (Nl < Noting the initial condition 
Qo(0) = 1 and Qo{L) = (i G N), wc find 



1 — (1 — a)z , , 



1 



M,(L + 1) 



(l-a)/3z ' ' (l-a)/3z' 
1 - (1 -a)(l -;3)z 



(1 - ol)^z 



-M,{L) 



-M,(L-1) 



(l-a)/3 
( ^■.XM,{L) + YM,{L-\)). 



(65) 



(66) 



The solution to the recurrence formula (l66l) is 



M,{L) = 



A 



A+- A_ 



A^ 



A, - A_ 



(M,(l) - A_M,(0)) 



(67) 



^^'^^ where A± = are the solutions to A^ 



^± ~ 2(l-a)/3z 

A A + Y with X and Y as defined in Eqn. Due 
to the condition |Afz(L)| < , the "initial condition" 
must be restricted as 

M,{l)-A^M,{0) = 0. (68) 

(Note that < |A_| < 1< |A+|.) Thus we find 



M,(L) = A4(0)A^ = 



l-{l-a){l + l3)z + r 



A^. (69) 



The generating function of the generating function is cal- 
culated as 



L>0 



M,(0) 



(70) 



and that of the average system length as 



t>o 



t>0 L>0 



-I + (1 + a - P + al3)z + r 



(71) 



4=1 2(1 -z)2 

The asymptotic behavior of (Lf) is determined by the 

1-I-/3' 



degree of the singularity z — 1 of 5^. When a < we 



find 

(1 - z)S, 
and (Lt) converges as 
(Lt) - 



1 /3 — a — a/3 

(t —i' oo). 



(72) 



(73) 



Again this limit value agrees with the stationary value 



© with p = 1. When a > we find 



(1 - z)^S, 
and {Lt) behaves as 



a- P + aP, 



{Lt) = 



a 



1+^ 



1 



t + o{t) 



t + o{t) (t->oo). 



When a = we find 



(l-^)*5. 



z-i-l 



1 + /?' 



(74) 

(75) 
(76) 

(77) 
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and {Lt) behaves as 



(i^.)-2^/lAl+,(v7) (t^oo). (78) 



V. CONCLUSION 

We have investigated the dynamical properties of the 
EQP, a queueing process with excluded- volume effect. 
The model can be interpreted as a TASEP with varying 
length. Using generating function techniques and a phe- 
nomenological domain wall theory we have derived ana- 
lytical predictions for the time-dependence of the number 
of particles {Nt) and the average system length {Lt). 

We found that the two phases observed previously can 
be divided in subphases. The convergent phase, where 
the system length remains finite, consists of high-density 
and maximal current subphases. The same is true for 
the diverging phase, where the system length becomes 
infinite in the long time limit. 

By comparing with Monte Carlo simulations it was 
found that the predications of the domain wall theory for 
the dynamical behavior are at least qualitatively correct, 
i.e. {Lt) and {Nt) converge or diverge linearly in time. 
Moreover they are in good agreement even quantitively 
in the convergent phase. In the divergent phase, the pre- 
dicted velocity for {Nt) appears to be correct whereas 
deviations from the domain wall theory can be observed 
for {Lt). 

For p — 1, we derived exact analytical results for the 
behaviors of {Nt) and {Lt) by using the generating func- 
tion method. We showed the linearity of their dynamics 
in the divergent phase as predicted by the domain wall 
theory. We found diffusive behavior on the critical line 
as well. 

The simple approach presented here does not provide 
a good expression for the velocity V of the system length 
{{Lt) ^ Vt) in the divergent phase. Here further analy- 
sis of the detailed density profiles in the divergent phase 
may be helpful. A first step would be the numerical de- 
termination of V to get a better understanding of its de- 
pendence on the parameters a, /3 and p. Another way to 
settle the problem is extending the exact result to p < I 
case, which may be difficult but very worthwhile. 



Appendix A: Stationary State for Inhomogeneous 
Injection Case 

Here we consider the stationary state for a general- 
ized model where the entry probability depends on the 
system length. A new particle enters the system with 
probability if the leftmost occupied site is L, or ao 
if there is no particle on the chain. The stationary state 
of this generalized model can be written in the following 
matrix product form with the same matrices and vectors 
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FIG. 8: Example for transitions from 1011 at time i — 4 to 
the states in which the system length is 5 and the 4th site is 
empty (left) or occupied (right) at time t. 



{D,E, {W\ and \V)): 

T-rL-l 

PiTL • • • ri) = - ,ir""' {W\Xr, ■■■Xr, \V), 



nL-l 



(Al) 

{W\D{D + E)''-'\V). 

(A2) 



This can be proved in the same way as in the homoge- 
neous case aL — a, see 

Appendix B: Derivation of Eqn. (1641) 



(Bl) 



To derive Eqn. ((64)) . we show 
PQtiL) 



Prob 



Prob 



the system length is L at time t 
A (L — l)th site is empty at time t 

{l~p)Qt{L) 

the system length is L at time t 
A (L — l)th site is occupied at time t 



(B2) 



for L e N>2 

In the deterministic hopping case p = \ with the initial 
condition 0, the sequence 00 (except the infinite number 
of O's left to the leftmost particle) never appears, and 
holes necessarily "hop" leftward. Thus the first site must 
be occupied by a particle at time < — L + 1 if the system 
length at time t is L. The local state (empty or occupied) 
of the {L — l)th site at time t depends only on whether 
the first particle exits or not at time t — L + 1. Let the 
number of particles be n at time t — L + 1, and label the 
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particles by natural numbers as in Fig. [8] We can show 
by induction that the £th particle is on the ^th site at 
time t — i + ^. Thus we find that L ~ n new particles 
should enter the system during [t — L + 1, t — 1]. 

We obtain (the summation of) the transition probabil- 
ity from a configuration 5*^*^^+^-' at time t — L + lto the 
states in which the system length is L and the (L — l)th 
site is empty or occupied at time t is given by 

t-i 



s(*),...,s(*-^+2)es t'=t-L+l 

|s(*)|=i„(sW)_c_i=o (B3) 

^^^^ a^-"(l-a)"-i/3, 



L — n 



E 



s(*),...,s<*-^+2)es t'=t-L+l 
s(')|=i,,(s('))j^_i=i 



1 \ L-n|■^ „\n-l. 



(B4) 



L — n 



a^-«(l-a)"-^(l-/3), 



respectively, where W^(S'(*') ^ S'(*'+i)), \S^*">\ and 
{S^*^)j denote the transition probability from S*'* ^ to 
5'^* the length of the state 5*^* ^ and the local state 

L~ 1 



of site j. The binomial ^ ^ ^ J gives the number of 
possibilities for when the new particles enter the system. 



These equations lead to Eqs. (IBip and (jB2l) . 

If the length is L at time t + 1, there are the following 
three possibilities at time t: 

(i) the length is L — 1 and a new particle enters (with 
probability a), 

(ii) the length is L, the (L — l)th site is occupied, and 
no particle enters (with probability 1 — a), 

(iii) the length is L + 1, the Lth site is empty, and no 
particle enters (with probability 1 — a). 



Then we achieve Eqn. ([Ml) for L e Z>2. For £ = 1, the 
case (ii) is replaced by 

the length is 1, the particle at the rightmost site 
does not leave, and no particle enters (with proba- 
bility (l-a)(l -/?)). 
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